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A non-Markovian stochastic Schrodinger equation for a 
quantum system coupled to an environment of harmonic os- 
cillators is presented. Its solutions, when averaged over the 
noise, reproduce the standard reduced density matrix without 
any approximation. We illustrate the power of this approach 
with several examples, including exponentially decaying mem- 
ory correlations and extreme non-Markovian periodic cases, 
where the 'environment' consists of only a single oscillator. 
The latter case shows the decay and revival of a 'Schrodinger 
cat' state. For strong coupling to a dissipative environment 
with memory, the asymptotic state can be reached in a finite 
time. Our description of open systems is compatible with dif- 
ferent positions of the 'Heisenberg cut' between system and 
environment. 



The dynamics of open quantum systems is a very timely 
problem, both to address fundamental questions (quantum 
decoherence, measurement problem) as well as to tackle the 
more practical problems of engineering the quantum devices 
necessary for the emerging fields of nanotechnology and quan- 
tum computing. So far, the true dynamics of open systems 
has almost always been simplified by the Markov approxima- 
tion: environmental correlation times are assumed negligibly 
short compared to the system's characteristic time scale. 

For the numerical solution of Markovian open systems, de- 
scribed by a master equation of Lindblad form 



-TT/'< = -' [// /'■] + ^ ^ ([L m pt, + PtL^]) 



(1) 



(where pt denotes the density matrix, H the system's Hamil- 
tonian and the operators L m describe the effect of the envi- 
ronment in the Markov approximation), a breakthrough was 
achieved through the discovery of stochastic unraveUings |l]J2j . 
These are stochastic Schrodinger equations for states 4>t[z), 
driven by a certain classical noise Zt with distribution func- 
tional P(z). Crucially, the ensemble mean M[. . .] over the 
noise recovers the density operator, 
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Hence, the solution of eq.(Q) is reduced from a problem in the 
matrix space of p to a much simpler Monte Carlo simulation 
of quantum trajectories ipt(z) in the state space. 

For the Markov master eq. ([[]), several such unraveUings are 
known. Some involve jumps at random times ji| , others have 
continuous, diffusive solutions Q]. They have been used ex- 
tensively over recent years, as they provide useful insight into 



the dynamics of continuously monitored (individual) quan- 
tum processes ||, or into the mechanism of decoherence Q. 
In addition, they provide an efficient tool for the numerical 
solution of the master eq.(|l]). It is thus desirable to extend 
the powerful concept of stochastic unraveUings to the more 
general case of non-Markovian evolution. 

The simplest unraveUings are linear stochastic Schrodinger 
equations. In the Markov case (|l|), for a single L, the linear 
equation 



d -ipt = —iHipt + Lipt o zt — \-L} Ltjjt, 



(3) 
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provides such an unravelling, where, z t is a complex-valued 
Wiener process of zero mean and correlations M[z^z s ] — 
S(t — s), M[z t z s ] = 0, and where o denotes the Stratonovich 
product 0. 

However, eq.(^) is of limited value, since the norm HV'tO 8 )!! 
of its solutions tends to with probability 1 and to infinity 
with probability 0, such that the mean square norm is con- 
stant. To be really useful, one should find unraveUings in 
terms of the normalized states 
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which requires a redefinition of the distribution of the noise 
P(z) -» P t {z) = \\ip t (z)\\ 2 P(z) H so that eq.(§) remains valid 
for the normalized solutions: 
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Now (^) can be interpreted as an unravelling of the mixed 
state pt into an ensemble of pure states. For the Markov 
unravelling (^|), the normalized states satisfy the non- linear 
Quantum State Diffusion (QSD) equation [^|: 



-i>t = -iHfa + {L- (L) t )$t o (zt + (L f ) t ) 



(6) 
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where (L) t = (ipt\L\ipt) ■ Contrary to eq.Q, eq.(^|) provides 
an efficient Monte- Carlo algorithm for the numerical solution 
of § 0§. 

In this Letter we present for the first time a nonlinear non- 
Markovian stochastic Schrodinger equation that unravels the 
dynamics of a system interacting with an arbitrary 'environ- 
ment' of a finite or infinite number of harmonic oscillators, 
without any approximation. In the Markov limit, this unrav- 
elling reduces to QSD (^) and will therefore be referred to as 
non-Markovian Quantum State Diffusion. Other authors have 
treated non-Markovian open systems effectively with Marko- 
vian unraveUings: either the system has in fact been influ- 
enced by a second Markovian environment in addition to the 
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original non-Markovian one or, alternatively, fictious modes 
have been added to the system Q. In our approach, the 
'system' remains unaltered, and the unravelling is genuinely 
non-Markovian. 

Below we summarize the general theory, which will be pre- 
sented in detail elsewhere |8|, and we present four examples: 
First, we consider a 'measurement-like' environment. Then, a 
dissipative environment with exponentially decaying environ- 
ment correlations is discussed. Remarkably, here the asymp- 
totic state can be reached in a finite time. In the third exam- 
ple we consider an 'environment' consisting of only a single os- 
cillator. This example is thus periodic, that is extremely non- 
Markovian. It shows the decay and revival of a 'Schrodinger 
cat' state. Finally, the fourth example shows that the de- 
scription of a subsystem in terms of non-Markovian QSD is 
independent of the 'Heisenberg cut', that is independent of 
where precisely the boundary between system and environ- 
ment is set. 

Our starting point is the non-Markovian generalization of 
the linear stochastic equation (til), derived in M, 
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ipt = — iHipt + LiptZt — / a(t,s)—^ds, 
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which unravels the exact reduced dynamics of a system cou- 
pled to an environment of harmonic oscillators. Here, zt is 
colored complex Gaussian noise of zero mean and correlations 
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a(t,s), M[z t z a ]=0. 
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The Hermitian a(t,s) = a*(s,t) is the environment correla- 
tion function |j |)| . The functional derivative under the mem- 
ory integral in (pp indicates that the evolution of the state ipt 
at time t is influenced by its dependence on the noise z a at 
earlier times s. In |§| we show that it amounts to applying an 
operator to the state, 

—iPt = 6(s,t,z)iJt, (9) 

where the explicit expression of 0(s,t,z) can be determined 
consistently from eq.(fij). 

Just as in the Markov limit, to be really useful, one has to 
find the corresponding non-linear non-Markovian QSD equa- 
tion for the normalized states (JiJ) . This quite elaborate deriva- 
tion can be found in Isf and leads to 
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iPt = -iH^ t + (L- {L) t )*p t zt (10) 



«(t,s) (AL t 6(s, t, z) - (AL t 6(s,t,5)) t ) ds^ t , 



which is the basic equation of non-Markovian QSD. Here, Zt is 
the shifted noise zt = Zt + J a a*(t, s)(L^) s ds, and for brevity 
we use AL 1 " = V - (L f ) t . 

Let's turn to concrete examples of non-Markovian QSD 
(^p|). First, we consider an environment modeling energy mea- 
surement: L = iJ — H. It is easy to prove that O — H in 
(|), and hence © reads 



■^ t = -iH4, t ~{H 2 - (H 2 ) t )$t a(t,s)ds 



(11) 



+ (H- (H) t )ipt I z t + J a(t, s)*{H) 3 ds + J a(t, s)ds{H) t 



Notice that indeed, flll| ) reduces to the Markov QSD equation 
(|) for a(t,s) -> S(t-s). 

If the correlation a(t, s) decreases fast enough, the asymp- 
totic solution of ([ll]) is an eigenstate <j> n of H, reached with 
the expected quantum probability |(0 n |-!/)o)| 2 . Numerical so- 
lutions of (0) for the 2-dimensional case H = ^a z and expo- 
nentially decaying correlation are shown in Fig. la (solid lines). 
The asymptotic state is either the 'up' or the 'down' state. 
The ensemble mean M[{a z )] remains constant (dashed line) 
as expected from the analytical solution (dot-dashed line). 
Note, however, that if the environment consists of a finite 
number of oscillators, represented by a quasi-periodic corre- 
lation function a(t,s), such a reduction to an eigenstate will 
not occur (see our third example). 

As a second example, we consider a dissipative spin with 
H — ^<7z, and L = A<r_. We choose exponentially decaying 
correlations a(t,s) = l e -~i\ t - s \- in ( t - 3 ) with an environmen- 
tal central frequency SI and memory time 7 _1 . The non- 
Markovian QSD equation (FlOI) reads pi 
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+ \(a--(o-)t)i>t\zt + \ J a(t,s)*{a+) s ds+ (a + ) t F(t) 

with F(t) determined from 
d 



dt 



F{t) = --yF(t) + i(u - Q)F(t) + \F(tf + ^ (13) 



and initial condition F(0) = 0. The equation for F(t) can 
be solved analytically (j. It is worth mentioning the case 
of exact resonance, ui — Q. Two regimes should be distin- 
guished. First, when 7 > 2 A 2 (short memory compared to 

coupling strength), F(t) tends to ^7 — \J 'i 2, — 27A 2 ^ /(2A). 

Hence, for large 7, one recovers Markov QSD (^|). For longer 
memory times or stronger coupling, 7 < 2A 2 , things are very 
different: F(t) diverges to infinity when the time t approaches 

t c = (jr + 2arctan(7/^2A 2 7 - 7 2 )) /^/2A 2 7 -7 2 . All re- 
alizations tpt(z) reach the down state in a finite time and 
remain there! In Fig. lb we show quantum trajectories from 
([[2]) (solid lines), their ensemble mean value M[{a z )t] (dashed 
line), and the analytical mean value (dot-dashed), which is 
almost indistinguishable. The reduction time in this case is 
Lot c = |7T ~ 4.71. This is the first example of a continu- 
ous quantum state diffusion that reaches its asymptotic state 
in a finite time, which was proven impossible for Markovian 
diffusions {[(J. 

Our third example is a harmonic oscillator coupled to a fi- 
nite or infinite number of oscillators initially in their ground 
states. Here, the non-Markovian QSD eq.(|10|) takes the same 
form (|l2]), where the Hamiltonian is uja^a and where <r_ (0 + ) 
has to be replaced by the annihilation (creation) operator a 
(aJ). The resulting equation preserves coherent states. More 
interesting is the case of an initial superposition of two sym- 
metric coherent states, known as a 'Schrodinger cat' If 
the environment correlation a(t,s) decays, so does the 'cat'. 
If, however, the environment consists of only a finite number 
of oscillators, then the 'cat' will first decay, due to the local- 
ization property of QSD, but since the entire system is quasi- 
periodic, the 'cat' will then revive! In Fig. 2 we show contour 
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plots of the evolution of the Q-function of such a 'cat', in the 
extreme case where the environment consists of a single os- 
cillator (a(t,s) — e -!n ''~ s '). Apart from an overall spiraling 
motion due to the 'system' Hamiltonian, the 'cat' state first 
decays but later revives. Our non-Markovian QSD equation 
thus provides a nice illustration of proposed experiments on 
reversible decoherence jL2j. 

As a last example we consider a case where the split be- 
tween system and environment can be shifted naturally be- 
tween two positions, see Fig. 3. A spin (Hamiltonian Hi) and 
a distinguished harmonic oscillator (H2) are linearly coupled 
(i/12). Moreover, the spin is coupled {Hi) to a heat bath 
(H env ) at zero temperature. We can either consider the quan- 
tum state of the spin-oscillator system coupled to a heat bath, 
or the quantum state of the spin coupled to a heat bath and 
coupled to the distinguished oscillator. In the first case, we 
can apply the Markov QSD description, ie a family of spin- 
oscillator states ipt(£) indexed by a complex Wiener process 
£t. In the second case, using non-Markovian QSD, we have a 
family of spin states <j> t (^,z) indexed by the same £t and also 
by the noise z t , due to the distinguished oscillator 'environ- 
ment' with correlation M[z*z s ] = e - iu ^t-s) _ 

Let us study a shift of the 'Heisenberg cut': compare the 
states 4>t{£,, z.) of the spin averaged over the noise z t with the 
mixed state obtained by tracing out (Tr2) the oscillator from 
the spin-oscillator states ipt(t;)- We prove in || that the states 
corresponding to both descriptions are equal: 



M z [\<f> t (£, *)> {<j> t (f , z) \] = TV 2 (| (£)> (f ) I )■ 



(14) 



This illustrates the general fact that non-Markovian QSD at- 
tributes stochastic pure states to a system in a way which 
depends on the position of the Heisenberg cut, but which is 
consistent for all possible choices of the cut. 

In conclusion, we present the first non-Markovian unrav- 
elling of the dynamics of a quantum system coupled to an 
environment of harmonic oscillators, which can thus be sim- 
ulated by classical complex noise. In the Markov limit, stan- 
dard Quantum State Diffusion is recovered. We emphasize 
that non-Markovian QSD (|lo| ) reproduces the true evolution 
of the system taking into account the exact unitary dynam- 
ics of system and environment The power of this new 
approach to open quantum systems is illustrated with four 
examples. For measurement-like interactions, reduction to 
eigenstates takes place whenever the environment correlation 
function decreases fast enough. For dissipative interaction 
with a heat bath at zero temperature, the ground state may 
be reached in a finite time. The third example is an appli- 
cation to the most extreme non-Markovian case: two linearly 
coupled oscillators, one of them playing the role of the 'envi- 
ronment'. We see the decay and revival of a 'Schrodinger cat' 
state. Finally, the last example illustrates that unravellings 
corresponding to different positions of the 'Heisenberg cut' 
between system and environment are mutually compatible. 
Most of these features are entirely new and have no coun- 
terpart in any Markov unravelling. Hence, non-Markovian 
unravellings represent a promising route to open systems, as 
for instance to quantum Brownian motion ^| . Moreover, our 
approach represents a new efficient tool for the numerical sim- 
ulation of quantum devices, whenever non-Markovian effects 
are relevant [filJli 
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FIGURE CAPTIONS 



FIG. 1. Non-Markovian quantum trajectories (solid 
lines) for a spin | system H — j<t z with an ex- 
ponentially decaying environment correlation a(t, s) — 
% exp (— 7|t — s\ — iQ(t — s)), where we choose 7 = w. The 
ensemble mean value over 10000 runs (dashed line) is in very 
good agreement with the analytical result (dot-dashed line). 
We show (a) a measurement-like interaction L — \a z with 
f2 = 0, and (b) a dissipative interaction L — Acr_ on resonance 
Q — ui, where each trajectory reaches the ground state in a 
finite time u)t c = §7r « 4.71. In both cases we choose a cou- 
pling strength A 2 = u> and an initial state |V>o) = 3| t) + 2| J.). 

FIG. 2. Reversible decoherence of an initial symmetric 
'Schrodinger cat' state \ipo) = \a) + \ — a) with a = 2. The 
contour plots show the Q-function of a non-Markovian quan- 
tum trajectory of a harmonic oscillator (u>) 'system', coupled 
to just a single 'environment' oscillator (fi = 0.5w). The cou- 
pling strength between the two oscillators is 0.1a;, and the 
time step between two successive plots is 0.47/w. 

FIG. 3. 'Spin - single oscillator - heat bath' system. First, 
we consider the 'spin - single oscillator' as the 'system' with 
state i/Jt(£), coupled to the heat bath with noise £ t . Alter- 
natively, we can consider the 'spin' as the 'system' t (£,z), 
coupled to the 'single oscillator + heat bath' environment 
(noises (£t,2t)). In non-Markovian QSD, both descriptions 
are possible and lead to the same reduced spin state. 
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